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Abstract 

We have found the exact ground state for two electronic models on a linear 
chain. The first model describes half-filling electron system in the ferromagnet- 
antiferromagnet transition point. In the singlet ground state the spin correlators 
show giant spiral magnetic ordering with the period of a spiral equals to the sys- 
tem size. The second electronic model describes the point where the ground state 
has giant spiral off-diagonal long-range order and, therefore, superconducting. We 
suggest the formation of the ground state with giant spiral order (ferromagnetic or 
off-diagonal) as a probable scenario of the subsequent destruction of the ferromag- 
netism and the superconductivity. 

1 Introduction 

In recent years there has been a growing interest in studying systems of strongly correlated 
electrons in relation with high- Tc superconductivity. Because of the difficulty in dealing 
with the many-body problems the exact results are rare. It is well-known that some ID 
electron models can be exactly solved by the Bethe ansatz. However, a lot of ID quantum 
systems do not obey the Yang-Baxter equation, and thus are non-integrable. Another 
approach leading to exact results consists in construction of the exact ground state wave 
function for some quantum systems. Recently considerable progress in this problem has 
been achieved by using the so-called matrix-product (MP) form of the ground state wave 
function. It allowed to find exact ground state for various ID spin models Its origin 

can be traced back to the S = 1 spin chain model 0]. For higher dimensional spin and 
electronic systems, there are also some methods for the construction of the exact ground 
state wave function 0-0. 

There is a class of the ID quantum spin models describing the ferromagnet - antifer- 
romagnet transition point, for which the exact ground state wave function has been found 
in 1^, |T0[. The singlet ground state wave function at this point has a special recurrent 



form, and for special values of model parameters it can be reduced to the MP form or to 
the RVB form ||^ . Spin correlations in the singlet ground state show giant spiral magnetic 
structure with the period of the spiral equals to the system size. On the antiferromagnet 
side of this point the ground state can be either gapless with the algebraic decay of spin 
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correlations or gapped with the exponential decay of correlations [0 . So, this model 
describes the boundary between the ferromagnetic phase and the singlet phase without 
long-range order. 

In this paper we present the singlet ground state wave function of this spin model in 
another form, which can be easily generalized for the ID electronic models. Then we con- 
sider two ID electronic models. The first model describes half-filling electron system at the 
point where the singlet and ferromagnetic states are degenerate. The exact calculation of 
the correlation functions in the singlet ground state shows the same giant spiral magnetic 
ordering as for the original spin model, while all other correlations vanish in the thermo- 
dynamic limit. The second electronic model describes the boundary on the phase diagram 
between the superconducting phase with off-diagonal long-range order (ODLRO) and the 
non-superconducting phase. The correlation functions in the ground state of this model 
show the giant spiral off-diagonal long-range order. We presume that in one-dimensional 
systems the destruction of the long-range order (ferromagnetic or off-diagonal) can be 
followed by the appearance of the ground state with the giant spiral order. 

We generalize this form of the wave function for the electronic ladder model. This 
model possesses both the giant spiral spin order and the giant spiral ODLRO in the ground 
state. Therefore, this electron ladder model describes the boundary points on the phase 
diagram between four different phases: two singlet phases with and without ODLRO, and 
two ferromagnetic phases with and without ODLRO. For some special cases the ground 
state wave function can be reduced to the usual MP form. 

The paper is organized as follows. In Section 2 we construct the exact singlet ground 
state for the quantum spin model. In Section 3 two electronic models with exact ground 
states are considered, and the correlation functions are exactly calculated. Section 4 gives 
a brief summary. In Appendix a technique for the calculation of correlators is developed. 



2 Quantum spin model 

First, we consider the s = | spin chain model with nearest- and next-nearest neighbor 
interactions given by the Hamiltonian 

N 11^ 1 

H = - ^(S, ■ S,+i - -) + - 5:(S. • - -) (1) 

i=l 1=1 

with periodic boundary conditions and even N. 

This model describes the ferromagnet-antiferromagnet transition point where the fer- 
romagnetic and singlet states are degenerate. The Hamiltonian (|l]) has been considered 
in [|], |T0[, where the singlet ground state wave function was constructed in two different 
forms. In this paper we represent another form of this singlet function, which allows us to 
generalize this function for the electronic model and to develop a technique to calculate 
correlators. 

The singlet ground state wave function for the Hamiltonian (p can be written as 
follows: 

^0 = ^0^, ^ = (Oblc/i® •..®^w|06) (2) 

where 

9^ = b^ |T). + Hi). (3) 
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Here we introduced one auxiliary Bose-particle 6"*", b (the Bose operators 6^, 6 do not 
act on spin states ||)^ and and the Bose vacuum |Ofe). Therefore, the direct product 
Qi® 92® ■■■ ® Qn is the superposition of all possible spin configurations multiplied on the 
corresponding Bose operators, like h'^hhh'^ ... |tiit •••)• -Po is a projector onto the singlet 



state. This operator can be written as 



2tt 



da 



2n 



dp / sin 76/7 e""" e 

JO 



(4) 



where S^^^^ are components of the total spin operator. 

This form of wave function resembles the MP form but with infinity matrix, which is 
represented by Bose operators. Therefore, we have to pick out the (Ofc|...|Ofe) element of 
the matrix product instead of usual Trace in the MP formalism because Trace is 

undefined in this case. The function contains components with all possible values of 
spin 5" (0 < S* < N/2) and, in fact, a fraction of the singlet is exponentially small at large 
A^. This component is filtered out by the operator Pq. 

In order to show that \E'o is the ground state wave function for the Hamiltonian (|I|), 
let us represent the Hamiltonian (|l]) as a sum of Hamiltonians hi of cells containing three 
sites 

N 

H = Y.h^, (5) 



where 



hi — — — (S^Sj-i 



i+2 



1,..N 



4' 



■4) -2(^^+1^-^+2-4) + 4' 

One can easily check that each cell Hamiltonian hi for i 
acting on gi O gi+i gi+2- 

Since each hi is non-negatively defined operator, then \E' is the exact ground state wave 
function of an open chain: 



2 gives zero when 



N-2 



H, 



open 



As it was mentioned above, the function \E' contains components with all possible values 
of total spin S, and, therefore, the ground state of the open chain is multiply degenerate. 
But it can be proven (as it was done in [|iy]), that for the cyclic chain (0) only singlet and 
ferromagnetic components of have zero energy. Therefore, for the cyclic chain (|l|) \E'o is 
the singlet ground state wave function degenerate with the ferromagnetic state. 

The exact calculation of the norm and spin correlation function (SjSj_|_/) (see Appendix) 
in the singlet ground state (^) results in the following expressions: 



(^ol^o) 



d 



N 



d^ 



N 



^cos^d) 



(^o|S.S,+,|^o) 



d 



N-l-2 



3cos(e-C)^ 



N-l-2 



COS 



(6) 



(7) 



It can be shown that in the thermodynamic limit equations 

{SiSi+i) = ^cos I'^j 



«=C=o 

j) and (1^) result in 
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So, we reproduce the result obtained in |T0[ that in the thermodynamic hmit a giant 
spiral spin structure is realized, with the period of spiral equal to the system size. 

3 Electronic models 

Now we will construct electronic models by generalization of the wave function (|^): 

^'o = Ps=o Trj (Ofel gi®g2®...®gN \%) (9) 

where 

9^ = h^\])^ + h\[), + cU^\2)^ + f\Q);) (10) 

with 6"*", h and / are the Bose and the Fermi operators respectively, |0f,) is the Bose 
vacuum, and c is a parameter of the model. Here we also denote an empty site by |0), a 
site occupied by one electron by |t) and ||) and a doubly occupied site by |2). So, the 
product gi ® ...gN is the operator in the Bose 6"^, h and the Fermi /"*", / spaces. We pick 
out the (Ofe|...|Ofe) element in the Bose space, which can be written as: 

(Ofcl gi®g2®-<^gN |0,) = 0o + (/+/ - //+) <Pi (11) 

and then we take the Trace over the Fermi operators: 

Trf{%\ gi®g2® - ® gN \0b) = 20o 

The projector Ps=o filters out the singlet component from the function 0o- 

Thus, we have a singlet wave function \l/o describing the state with one electron per 
site. 

In order to find the Hamiltonian for which the wave function (^) is the exact ground 
state wave function, let us consider what states are present on the two nearest sites in 
the wave function (P). One can easily check that there are only 9 states from the total 16 
states in the product g^ (g) gfj+i. They are 

ITT), lii), iTi + iT), 120 - 02), III - It) 120 + 02), 
110 -0 1), IT 2-2 1), liO-Ot), 112-2 I) (12) 

The elementary Hamiltonian for which all these states are the exact ground states 
can be written as the sum of the projectors onto the 7 missing states \ipk) with arbitrary 
positive coefficients 

7 
k=l 

For c > 1 the most simple form of this Hamiltonian is: 

N 

H = Y.K^+l (13) 
i=l 

3 4 

hi,2 = l-4:Si-S2 + 4:{l--)rilril + -rii-ri2 

+ H(ct<7C2,<7 + C^<xCl,<7)(l - -a - ^2 -a) 
a 

2 

+ — J2('^taC2,a + C\„CX,„) (^i _^ - n2-af 
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Here c^^., Ci^a are the Fermi operators, ni^„ = cf^^Ci^a, and the SU{2) spin operators are 
given by Sf = c^^ Ci^i, = c^^ and = |(?^i,t — We also use here r] operators: 

Vl = 4,1 <T' Vi = Ci,^ Ci^h Vi = ^ , 



which form another SU{2) algebra ]T5|, |l6l, and rji ■ 772 is a scalar product of pseudo-spins 
rji and 772- 

The Hamiltonian (|13D does not conserve the total number of empty and doubly occu- 
pied sites because of the last term in the elementary Hamiltonians , in contrast to 



the models considered in 16 



Each elementary Hamiltonian {i = 1, ..A^ — 1) acting on functions 00, 0i gives 

zero, since all the states (|1^) are the eigenstates of /ij^j+i with zero energy, while the 
energies of all other states at c > 1 are positive. Therefore, the functions 0o and 0i are 
the ground state wave functions of the open chain: 



N-l 



i=l 



To determine the degeneracy of the model (p!4D , we need to classify the functions 0„. 
Analogously to the spin model (^, the functions 0o cind 0i contain components with all 
possible values of total spin S. Therefore, 0o contains multiplets with S = 0, ...N/2 and 0i 
contains components with values of the total spin S = 0, ...N/2 — 1 (0i does not contain 
the ferromagnetic component, since at least two sites in 0i are non- magnetic |0) and |2)). 
So, -|- 1 multiplets are degenerated for the open chain. 



But for the cyclic model (13) it can be proved that only three multiplets are the ground 
states: singlet state with the momentum p = (singlet component of (po), the trivial 
ferromagnetic state with p = tt, and the state with 5* = N/2 — 1 and p = tc (which is the 
component of (pi with S = N/2 — 1). The last state with 5* = N/2 — 1 can be written as: 

N 

^7V/2-i = E«TC.,i + c+T^^i) n <i |0) 

i<j n=l 

Thus, the ground states of the electronic model (|13D with one electron per site are the 
singlet state, the ferromagnetic state, and the state with S = S'max — 1- 

It is interesting to note that the singlet wave function (|^) can be also written in the 
form: 

N 

^o = j:[^,j][k,l][m,n]... n<iio) 

n=l 

where 

[hj] = - 5+ + (c+^Cj,| - C+-fQ,|), 

and the summation is made for any combination of sites under the condition that i < j, 
k < I, m < n .... This form of the wave function is analogous to the RVB form found in 
[y] for the spin model (p. 

The norm and the correlators of the electronic model ([T3|) in the singlet ground state 
are calculated in the same way as for the spin model (Appendix). Therefore, here we give 
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the final results: 



(^0 

(^o|S.S,+, 



+1 



N 



COS 



21 



5=0 



I ptN-l-2 



& d 



Ql qN-1-2 
I f)/-N-l-2 



3 cos(^ - C)(l + cosh(c2^ + c\)y 



COS 



=C=o 



— c 



(cos(0 + cos(C))(cosh(c2^) + cosh(c\))^ 



2cos4(i±i: 



cosh(c^^ - c^C) ' 
2cos2(^) 



«=C=o 



d 



'N-2 



lN-2 



d^'' ^ \cos 



^1) 



(15) 
(16) 

(17) 

5=C=o 

(18) 
(19) 



5=0 



As it can be seen from Eg . (|T9D , the expectation value {r]ir]^_i), which determines the 
off-diagonal long-range order (ODLRO) |]T3|, does not depend on the distance I. But in 
the thermodynamic limit ODLRO vanishes: 



{eta Ci+l,a) 



0{ 



(SiSj+i) = - COS 



27r/ 



(20) 



though for finite systems all correlators (|T^-(|T^ are non zero. 

The second electronic model can be obtained by simply interchanging of the Bose and 
the Fermi operators in (PH]). So, the wave function of this model has the form: 



with 



^0 = Pv=o Trj (Ofel gi®g2® ■■■ ® Qn \0b) 



9^ = cin]), + f\^),) + b+\2), + b\0)^ 



(21) 



(22) 



The projector Pr,=o filters out the state with total 1] = = ^- Therefore, the function 
^'o has = 0, but it is not an eigenfunction of 5*^. Instead, it is an eigenfunction of r]'^ 
with rj = 0. 

The wave function (^If) can be also written in the form analogous to the RVB one: 

^o = T.[hj][k,l][m,n]...\0) 

where 

and the summation is also done over any combinations of sites under the condition that 
i<j,k<l,m<n .... 

Considering the product gi (g) gi+i one can find that there are only following 9 states 
on two nearest sites in the wave function (^Tj) 

|22), |00), 120 + 02), in -it), |20 - 02) + c2 Iti + it), 

IT + T), IT2 + 2T), liO + OT), Ii2 + 2i) (23) 
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The Hamiltonian for this model has the form which is similar to the previous one given 
by Eq.®: 

N 

H = Y.K^+l (24) 

i=l 

/ii,2 = l-4r7i-r72 + 4(l-4)5r52 + 4si-S2 

~ Y^{'^l,aC2,a + C^<^Ci,ct)(1 - _^ - ^2 -a) 

(7 

2 

+ — XI cr(c^;^C2,CT + C^^Ci,^) (rii _^ - ^2 -a)^ 

^ u 

This Hamiltonian for c > 1 is also a non-negatively defined operator and \E'o is the 
exact ground state wave function with zero energy. This Hamiltonian commutes with rf, 
but does not commute with S*^. Therefore, the eigenfunctions of the Hamiltonian (|^) can 
be described by quantum numbers rj and r^^. Making up the same analysis as for previous 
model we find that for the cyclic model (p^) the states with three different values of 77 have 
zero energy (as it was for model (0)). They are: one state with 77 = and momentum 
p = n ([21|), all states with i] = N/2 and p = 0: 



^N/2,r = iv^f^''^''\^) (25) 
and the states with 77 = N/2 — 1 and p = 0: 

^ N/2-1, = iv^r^'-''^' E«T4i - <Ad I 0) (26) 

i<j 

Therefore, for the case of one electron per site [i]^ = 0) the ground state of the model 
(|4D is three- fold degenerate. 

The correlation functions in the ground states (^) and (|26|) obviously coincide with 
each other in the thermodynamic limit and for half filling case [rj^ = 0) they are: 

«,c,+,,.) = 0(l) (SA+,) = 0(i^) 
VM+^ = O(^) (v-vti) = \ + O(^) (27) 



The existence of ODLRO immediately follows from the form of the wave functions ( [25| ) 
and (p^). 

The correlation functions in the ground state ( pi] ) have similar forms as in Eqs. ([TB|)- 
(p!9D, and in the thermodynamic limit they reduce to 

«,Q+,,.) = 0(1), (S,S,,^,) = O(^), 



(vivU) = 2 (vtvUi) = I COS (^) (2J 



6 \N ^ 

The giant spiral ordering in the last equation implies the existence of ODLRO and, 
therefore, the superconductivity [|T^ in the ground state (pTf). We note that though all 
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Figure 1: The two- leg ladder 



three ground states of the model ( p^ are superconducting, the properties of these wave 
functions are essentially different. Let us consider density - density correlator {niUi+i). For 
the wave functions ( ^5l) and (26) in the thermodynamic limit this correlator decouples: 
{riirii^i) = {rii) (rii^i) = 1. But for the wave function (PD it equals to {riiTii^i) = 1 + 



■ cos 



It is interesting to note that another model having the ground state wave function pl| ) 
with c = and the same spiral ODLRO (^) can be obtained from the model (|lD by simple 
replacing of operators S by rj: 



N 



N 



(29) 



1=1 



=1 



The direct analogy of this model to the spin model (|l]) results in the conclusion that the 
model ( p9|) describes the boundary point on the phase diagram between superconducting 
and non - superconducting phases where the off-diagonal long-range order is destroyed. 
We suppose that the model ( ^4|) also describes such a point. Thus, the wave functions 
(^ and (|21| ) are the ground states for the ID electronic systems in the boundary points 
between the phases with and without long-range order (ferromagnetic for (|13D and off- 
diagonal for (P^)). And we suggest the formation of the ground state with the long-range 
spiral order like (1^) and (pSf ) as a probable scenario of the subsequent destruction of the 
ferromagnetism and the superconductivity. 

The proposed form of wave function can be further generalized for the electronic ladder 
model (Fig.l). The wave function for the cyclic ladder model containing 2N sites has the 
form: 

^0 = Ps=v=o {0\9i®92® ...®9n\0) (30) 
where each Qi corresponds to i-th rung of the ladder: 

9^ = ci (a+(2x - a+a) ITT), - a lii)^ + (a+a -x) in + it),) +C2 in -it), 

+ C3 (b+{2y - b+b) |22), - b |00), + {b+b - y) |20 + 02).) + C4 |20 - 02) . (31) 



with a"*", a and 6^, b are the Bose operators, |0) in (^) is the Bose vacuum of and fo"*" 
particles, and c,, x, y are the parameters of the model. P^^s=o is the projector onto the 
state with S* = r/ = 0. 

Let us first consider the case 03 = 04 = 0. In this case the wave function ( ^OD describes 
the spin ladder model depending on the two parameters C2/C1 and x. It can be shown 
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that this model coincides with that considered in |]T2| . It has the singlet ground state ( pO| ) 
degenerate with ferromagnetic state. The spin correlators in the singlet ground state show 
double-spiral ordering with small shift angle Ayj = between the two giant spirals 

formed on two legs of the ladder: 

(S„S„+2m) = -cos(^^ + (-irA<^j (32) 
For the cases of integer or half- integer x = j, which correspond to the special cases of 





(2ttI 


J cos 










(2tiI 


^cos 









the model [12], one can easily recognize in Eq.(RTf) the Maleev's boson representation of 



spin S = j operators: 

= a^{2j — a^a), S~ = a, = a^a — j 

Therefore, in these special cases the infinite matrices formed by the Bose operators a"*", a 
can be broken off to the size n = 2j + 1 and the wave function ( PPD is reduced to the usual 
MP form. The spin correlators in the special cases have the exponential decay. 

Now let us come back to the general case of the electronic ladder model (^). In 
order to find the Hamiltonian for which (|30D is the exact ground state wave function, one 
should consider what states are present in the \E'o on the two nearest rungs of the ladder. 
There are only 26 states from the total 256 states in the product Qi ® gi+i- Therefore, the 
Hamiltonian of the ladder model can be written as the sum of the projectors onto the 230 
missing states \ipk) with arbitrary positive coefficients A^: 

N 230 

H = J2^i,i+^^ hi,i+i = ^><k\y^k){fk\ (33) 

i=l k=l 

Unfortunately, we can not give the explicit form like ( [TB| ) and (^^ for this Hamiltonian, 
because it has a very cumbersome form. But we are able to determine some properties of 
the Hamiltonian (|33|). This Hamiltonian commutes with both S'^ and t]'^. It has multiply 



degenerated ground state: the state with S = rj = (|30| ) and all the states with S + ri = N 
have zero energy. Hence the electronic ladder model (|33D describes the boundary point 
between the phases with and without ferromagnetic and off-diagonal long-range order. 
The correlation functions in the ground state (|30| ) can be calculated with the use of the 
technique developed in Appendix. For the case Ci > C3 there is the same double-spiral spin 
ordering (^) as for the spin ladder model, while all other correlations are exponentially 
small. For the case ci < C3 the double-spiral ODLRO is realized. In the most interesting 
symmetric case Ci = C3, C2 = C4, x = y the system possesses both the giant spiral spin 
order and the giant spiral ODLRO: 

(S„S„+2/) = {r]nVn+2l) = ^ COS 



(S„S„+2«+i) = {rinVn+2i+i) = I COS + Aip] (34) 



I 

~N 



where Ao? = 
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Therefore, in this case the wave function ( pO] ) describes the boundary points on the 
phase diagram between the four different phases: the singlet phases with and without 
ODLRO, and the ferromagnetic phases with and without ODLRO. In the special cases 
when X (y) is integer or half-integer the spin (off-diagonal) correlations exponentially decay 
and the wave function in the corresponding Bose space can be represented in the MP form 
with finite matrices of the size n = 2x + 1 or n = 2y + 1. When both x and y are integers 
or half-integers the wave function (^) can be written in the usual MP form with the size 
of matrices n = {2x + l){2y + 1). 



4 Summary 

We have found another form of the singlet ground state wave function for the quantum 
spin model considered previously in [|, |10|. The special technique was developed for the 



exact calculation of the norm and the correlation functions. This form of the wave function 
allowed us to generalize it for two ID electronic models. 

The first model describes the half-filling electronic system in the ferromagnet - antifer- 
romagnet transition point when the singlet and the ferromagnetic states are degenerate. 
In the singlet ground state the spin correlators show the giant spiral magnetic ordering 
with the period of the spiral equals to the system size, while all other correlations vanish 
in the thermodynamic limit. 

The second electronic model in the half filling case has three-fold degenerate ground 
state. All ground states have off-diagonal long-range order and, therefore, are supercon- 
ducting. The calculation of the correlation functions shows that one of the ground states 
has a giant spiral ODLRO. 

The comparison of these electronic models with original spin model ITT], O leads us to 



the conclusion that these two electronic models describe the boundary points on the phase 
diagram between the phases with and without long-range order (ferromagnetic for the first 
and off-diagonal for the second model). Therefore, we presume that if the Hamiltonian 
of the ID quantum system commutes with operators forming the SU{2) algebra (it can 
be the spin S or the pseudo-spin t] operators), then the appearance of the ground state 
with the giant spiral order predicts the following destruction of the ferromagnetism or the 
sup er conduct i vity. 

We have briefly considered the generalization of the proposed form of the wave function 
for the electronic ladder model. The general case of this model has much more rich phase 
diagram than two first models. In some particular cases this model describes the boundary 
points on the phase diagram between four different phases: with and without ferromagnetic 
and off-diagonal long-range order. There are also some special cases of the electronic ladder 
model when the ground state wave function is reduced to the usual MP form. Besides, the 
proposed form of the wave function can be also generalized for the 2D case and different 
types of lattices. 
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Appendix 



Let us calculate the norm of the wave function \I'o (0) 

(^o|*o) = (^|Po|^) (A.l) 
Since the function \l/ has = 0, then the projector Pq in the Eq.(Q) takes the form [|T^ 

Po = \ r sin e'^^'e'"'^^ (A.2) 



2 







where z = tan ^, 2;' = sin ^ cos ^ and S'^^ ^ are the operators of the total spin. 
Therefore, the norm takes the form: 

1 N 



{^o\^o) = 77 rsmjdj{0,,Ot,\f[{gte''''e''''''g,)\0,,0,) 
2 Jo 

where 

g+^.s-^.'stg^ = (a+(T.| + aa.|)e-^re-'5,+ (6+|t), + 6|i)J 
= a'^b^ + (1 — z'z)ab + izab'^ + iz'a^b 

and a"*", a are the Bose operators. 

So, the norm can be rewritten as follows 

(v|/o|v|/o) = ^ rsin7d7(0|G^|0) (A.3) 

Z JO 

with |0) = |0a, Ob) is the Bose vacuum of and 6^ particles and 

G = u{a'^b'^ + ah) + if (a&+ + a+6) 

where u = cos 2 ? = sin |. 

Let us introduce the auxiliary function P{^): 

P{0 = (0|e««|0) (A.4) 

then, 

(0|G^|0) 



In order to find P{C,) we perform the following manipulations. First, we take the 
derivative of P{^): 

rIP 

= (0|Ge«^|0) = (0|e«^G|0) = u (0|a6e«^|0) = u (0|e«^a+6+|0) (A.5) 

Now we need to carry a^b^ over G in the last expression 

e«^a+6+ = e«^a+e-«^e«^6+e-«^e«^ (A.6) 
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It can be made by using the following equations: 

a"^ cos^ + {ub + ivb~^) sin^ 



e^'^a e~^^ = a cos ^ — {uh^ + ivh) sin ^ 



e^'^h^e ^'^ = 6^ cos^ + (wa + if a"*") sin^ 
e^^he'^^ = b cos ^ — {ua^ + iva) sin ^ 

Substituting Eqs.(|X:^), (^) into Eq.(^), we find 

(0|e«^a+6+|0) = usin^cos^ (0|e«^|0) + sin^ ^ {0\abe^^\0) 

The last equation can be rewritten as the differential equation on the P{^) 

—— = u sin4 cos 4 -T (?) + M sm 4 — — 



with boundary condition P(0) 
The solution of Eq. (|A.8| ) is 



P{0 



1 — sin 



Integrating Eq. (|A.3|) over 7, we obtain 

(^ol^o) = 
So, finally we arrive at 



1 r , d^P 

2 I ^^"^^^ 



d 



N 



cos^d) 



(^o|^o)-2^ tan| 



di^ 

4 (2^+^ - 1 
iV + 2 



JV+2I 



^+ gi.'s, 2S,^^, = M(a+6+ - a6) + if (a6+ - 



Here i?Ar are the Bernoulli numbers. 

To calculate the spin correlators we need to introduce operators: 

G_ = g+e"^^e''''^^S-gi = uab+ + iva+b+ 
Then, the correlator (SiS^+i) will be defined by 



g+ ^izs^ ^iz's^ g+g, ^ ^^+^ _^ -^^j^ 



(M/o|SiSi+i|M/o) = - / sin 7^7 (0| -G^G^G^G 

Z Jo 4 z 

(since (0|G'_...|0) = 0). 

The expectation values in Eq. (|A.12|) can be represented as 



(0| G,G^G,G 



N-l-2 



(0| 



|0) 
|0) 



al aN-l-2 



«=C=o 



5=C=o 



(A.7) 



(Ai 



(A.9) 



(A.IO) 



(A.ll) 



|0) (A.12) 



After the procedure similar to that for the norm and the integration over 7. 

gl QN^l-2 I 3cos(e-C)' 



(^o|SiSz+i|^o) 



9^' dC 



N-l-2 



COS 



(A.13) 
we obtain 
(A. 14) 



S=C=o 
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